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ABSTRACT 

In this paper we investigate the growth rates of Gortler vortices m a compressible flow 
in the inviscid limit of large Gortler number. Numerical solutions are obtained for 0(1) 
wavenumbers. The further limits of (i) large Mach number and (ii) large wavenumber with 
0(1) Mach number are considered. We show that two different types of disturbance modes can 
appear in this problem. The first is a wall layer mode, so named as it has its eigenfunctions 
trapped in a thin layer near the wall. The other mode we investigate is confined to a thin 
layer away from the wall and termed a trapped layer mode for large wavenumbers and an 
adjustment layer mode for large Mach numbers, since then this mode has its eigenfunctions 
concentrated in the temperature adjustment layer. We are able to investigate the near 
crossing of the modes which occurs in each of the limits mentioned. 


iThis research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No NAS1-18605 while the authors were in residence at the Institute for Computer Applications m 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1. Introduction 


Our aim in this paper is to investigate the growth rates of Gortler vortices for a com- 
pressible inviscid flow over an infinite cylinder in the limits of high Mach number and high 
wavenumber. This investigation is motivated by recent interest in the development of hyper 
sonic aircraft which might well be capable of reaching speeds in the order of 20 - 25 Mach. 
Real gas effects will certainly come into play at these speeds but for simplicity s sake we have 
not taken them into account in this paper. We have also assumed that Chapman s viscosity 
law holds for this fluid. 

The most obvious difference between Gortler vortices in incompressible and hypersonic 
flows is that the presence of a temperature adjustment layer, where the temperature decays 
rapidly to its free stream value, at the edge of the boundary layer enables hypersonic Gortler 
vortices to be concentrated well away from the wall. In the incompressible case, we know, 
from the work of Hall (1982a, b, 1983), and Denier, Hall, and Seddougui (1991), that unstable 
Gortler vortices are not localized within the boundary layer for order one Gortler numbers. 
For higher Gortler numbers the most dangerous Gortler vortices have wavelengths small 
compared to the boundary layer thickness and are trapped near the wall. Not surprisingly 
this situation does not change significantly for order one Mach numbers and this nonparallel 
problem has been discussed by Wadey (1990) and Spall and Malik (1989). In the latter 
two papers, the nonparallel equations were solved numerically following the approach of Hall 
(1983) and the main result obtained was that the growth rate of a Gortler vortex is a function 
of its upstream history. The numerical calculations of Wadey (1990) do suggest that as the 
Mach number increases the unstable Gortler vortices locate themselves towards the edge of 
the boundary layer. This view is supported by Hall and Fu (1989) whose main result was that 
the logarithmically small temperature adjustment layer at the edge of a hypersonic boundary 
layer can support Gortler vortices and the most dangerous wavelengths of the vortices are 
comparable with the thickness of this layer. Fu, Hall, and Blackaby (1990) have considered 
the influence of real gas effects and Sutherland's viscosity law on the Gortler instability in 
hypersonic flows. 

This present paper is restricted to linear regimes of vortex growth, for a detailed account 
of nonlinear regimes the reader is referred to the review article by Hall (1990). It is also 
worth noting that hypersonic boundary layers are susceptible to instabilities other than 
those caused by streamline curvature, such as Tollmien-Schlichting wave instabilities which 
have been discussed by Cowley and Hall (1990) and Smith and Brown (1990). Clearly any 
nonlinear investigation of Gortler vortices at hypersonic speeds must allow for the possible 
interaction of the vortices and these other types of instabilities. 

The present paper is concerned with the inviscid limit of Gortler vortices in a compressible 
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boundary layer. We find that two distinct modes exist and we are able to describe them 
in the separate limits of large free-stream Mach number and large disturbance wavenumber 
(although for the case of large free-stream Mach number the disturbance wavenumber will 
also be large). The first mode also exists in an incompressible boundary layer and is termed 
a wall layer mode here since its eigenfunctions are concentrated in a thin layer near to the 
wall. For the incompressible case, in the limit of large wavenumber, this mode is described by 
Denier, Hall, and Seddougui (1991). The second mode has its eigenfunctions concentrated in 
a thin layer away from the wall and therefore is referred to as a trapped layer mode in the large 
wavenumber limit. In the hypersonic limit this layer is precisely the temperature adjustment 
layer, a logarithmic layer where the basic temperature changes rapidly from its 0(M^ o ) value 
close to the wall to its 0(1) value at the edge of the boundary layer. Consequently, in the 
hypersonic limit this mode is termed an adjustment layer mode. 

The numerical solution of the equations governing the stability of Gortler vortices in the 
inviscid limit, discussed in Section 2, show that the growth rates of the infinity of solutions 
of the two modes described above, as functions of the disturbance wavenumber, appear 
to intersect, A similar near-crossing of modes is evident in the numerical results of Mack 
(1987) in the form of kinks in the neutral curves of the generalized inflection point modes 
for compressible flow over a flat plate. Asymptotic solutions for the infinity of solutions of 
the compressible Rayleigh equation, termed acoustic modes, in the hypersonic limit have 
been given by Cowley and Hall (1990), while the so-called, single, vorticity mode for large 
values of the Mach number has been discussed by Smith and Brown (1990). Cowley and 
Hall (1990) postulated that the near-crossing of the neutral curves could be described by 
a WKB description of the acoustic modes and the vorticity mode. This is precisely the 
method employed by Smith and Brown (1990) in their investigation of inviscid modes of 
instability for large Mach number flows. The results of Smith and Brown (1990) show that 
the discontinuous vorticity mode becomes continuous in the limit of large Mach number. 
Additionally, they ascertain that for large Mach numbers the acoustic and vorticity modes 
are separated by an exponentially small amount as was proposed by Cowley and Hall (1990). 

The objective of the present paper is to describe the wall layer modes and trapped layer 
modes present for inviscid Gortler vortices and to investigate their near-crossing which exists 
as outlined above. To this aim we follow the ideas of Cowley and Hall (1990) and Smith and 
Brown (1990) and consider a WKB analysis of the modes of interest. 

The layout of this paper is as follows. In Section 2 we derive the equation governing the 
structure of inviscid Gortler vortices in a compressible boundary layer and then discuss some 
numerical results of this equation. In Section 3 we take a closer look at the wall layer modes, 
those with their eigenfunctions trapped near the wall, in the hypersonic limit. We follow 
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this by looking at the adjustment layer modes, those whose eigenfunctions are concentrated 
in the temperature adjustment layer, in the hypersonic limit in Section 4. Then in Section 
5 we look at the wall layer modes for the high wavenumber limit.. In Section 6 we consider 
the trapped layer modes for the high wavenumber limit, these are the equivalent of the 
adjustment layer modes for the hypersonic limit. We then consider the near crossing of the 
different types of modes for both limits in Section 7 before finally summarizing our results 
in Section 8. 


2. Formulation 


Our aim in this section is to obtain the equation which determines the structure of a 
Gortler vortex in a compressible boundary layer. The boundary layer considered is that of 
a flow over the cylinder y* = 0, — oo < z* < oo so that the z -axis is a generator of the 
cylinder and y* measures the distance normal to the surface. The x -coordinate measures 
distance along the curved surface, which is supposed to have variable curvature (l/m)x(^ //) 
where m and l are length scales. The Reynolds number R, Gortler number G and curvature 
parameter 8* are defined by 


R = 


U^l 


G = 2RH\ 
l 


(2.1a) 

(2.16) 

(2.1c) 


m 


where Uoo is a typical flow velocity in the streamwise direction and v is the kinematic viscosity 
of the fluid. The Reynolds number is assumed to be large, whilst 8* is sufficiently small so 
that as 8* — i 0 the parameter G is fixed and of order one. We take the basic two-dimensional 
boundary layer to be of the form 


u = Uoo[u(X, T), R 2 v(X, T), 0] [1 -1- 0(R 2 )], 


( 2 . 2 ) 


where 


X- y-Ri 

X ~T' Y ~ l 


and T] is given by 


V = Y */X>, 


(2.3) 


We chose to look at a Blasius boundary layer by putting u = /'(??) where / satisfies 

2/"' + //" = 0, /(0) = /'( 0) = 0, /'(oo) = 1, (2.4) 


(2.5) 
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where Y* is the Howarth-Dorodnitsyn variable 

r ’ = £ r iY ’ (2 - 6) 

and we have taken a model fluid (i.e., the Prandtl number equal to one and Chapman’s vis- 
cosity law with C equal to one). For the case of a thermally insulating wall, the temperature 
is given by 

T = 1 + 5(7 - 1)A£(1 - P ), (2.7) 

where M, is the Mach number of the free stream and 7 is the ratio of the specific heats and 
will be taken to be 1.4 when it is needed numerically. Defining Z by 


Z = 


z'Rh 

I - ’ 


(2.8) 


we perturb ( 2 . 2 ) by writing 


u = U 00 [u + 6U(X,Y)E 1 , vR ~2 +8R~sV(X,Y)E 1 , <5iHw(X, y)£ a ][l + 0(iT*)], (2.9) 


where E\ = exp(zaZ) and { < 1 (see Hall (1982a) for further discussion of the above 
scalings). We similarly perturb the basic pressure, p, by putting 

p = p + 6R~ 1 P{X,Y)E 1 , (2.10) 

and the basic temperature by 

T = T 00 (T + ST{X,Y)E 1 ), ( 2 . 11 ) 


where is the free stream temperature. If we then introduce a growth rate £ and the 


scalings 

U{X,Y) = U(Y)E 2 , ( 2 . 12 a) 

V(X,Y) = G*V(Y)E 2 , ( 2.12 b) 

W(X, Y) = G 1 >W(Y)E 2 , (2.12c) 

T(X, Y) = T(Y)E 2 , ( 2 . 12 d) 

P(X,Y) = GP(Y)E 2 , (2.12c) 


where 


E 2 =exp(y Gif)(X)rlX ] j , 


( 2 . 12 /) 


and insert (2.9) - (2.11) into the governing equations and take the inviscid limit G->oowe 
get 


P0U . d (pV\ ia P w fipuT 
T dY\T ) T T 2 


(2.13a) 


4 


rs 

ufrU + V-^ = 0 , 


£ 

T 


u/?F + KuU 


1 1 u 2 K T 

2 T 2 


dP_ 

' dY' 


= -iaP, 

T 

dT 

up T + V qy=^ 

where K is a measure of the curvature. Taking p = 1 we find from (2.13) that 


V 


Ka 2 d (u 1 ' 

2 JudY\f t 


pa 2 u 


d ( fi du ' 
dY [TOY, 


+ $u 


d_ ( }_dV_" 
&Y { Toy, 


= 0. 


(2.136) 

(2.13c) 

(2.13d) 

(2.13e) 

(2.14) 


If we now change variables to rj, take for simplicity’s sake a value of X = 1 and put 


P - K*p, 


(2.15) 


in order to eliminate K from the equation we get after some rearranging 

r , 2 Tf* 


V" + V' 


2 T 

~¥ 


+ v 


' 2 T‘ - + -=- 


f T f‘ 


l) + *( T r_£ V 

') P 2 \ f 2). 


0, 


(2.16) 


with the boundary conditions 

V(0) = 0, (2.17a) 

V ~ e~ ar> as j] — > oo, (2.176) 

V‘ ~ —ae~ aT] as f] oo. (2.17c) 

Since (2.16) corresponds to an inviscid limit we cannot satisfy the viscous boundary condition 
V' = 0 at 7 / = 0. Equation (2.16) is the compressible generalization of equation (5.8) in 
Denier, Hall, and Seddougui (1991). Similarly to Denier, Hall, and Seddougui (1991) it can 
be seen that (2.16) and (2.17) have the exact solution 

V = f'e- a \ /? 2 = ^, (2.18) 

which is valid for all a and M <». However, in order to get other solutions it is necessary to 
solve this eigenvalue problem numerically. 

We solved (2.16) and (2.17) for a variety of Mach numbers. Figure 1 shows the first 
eleven modes for a Mach number of 2. These modes are not markedly different in shape 
from those of the form /3 2 = constant x a obtained in the incompressible case (again, see 
Denier, Hall, and Seddougui (1991)). Figure 2 shows solutions for a Mach number of 3.25 
and we can see that there are now significant differences compared with the incompressible 
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case. However, all of the modes continue to rise as a — + oo and are of the wall layer type, 
so-called because the eigenfunction is concentrated in a layer near the wall. In Figure 3 
where the Mach number is now 5, we can clearly see the appearance of two adjustment layer 
modes which tend to a constant value of (3 as a — *■ oo. We call these modes adjustment 
layer modes because their eigenfunctions are concentrated in the temperature adjustment 
layer in the hypersonic limit. (They will be referred to as trapped layer modes for 0(1) Mach 
numbers.) Each of the modes in this figure (apart from the exact solution) contributes a 
portion to the adjustment layer modes as they pass through the relevant positions and in 
doing so they change up to the next higher wall layer mode. If we follow the fourth mode 
on Figure 3 as a increases we can look at its eigenfunction and see how it changes. Initially 
when a = 0.2 we can see in Figure 4a that the eigenfunction is a widened version of the 
eigenfunction of the third wall layer mode with an additional lower peak in the temperature 
adjustment layer. In Figure 4b with a = 1.1 the mode is now in its flat section and the 
eigenfunction is indeed that of an adjustment layer mode. By the time a = 1.25, Figure 4c 
shows that the peak in the temperature adjustment layer has split into two. In Figure 4d we 
have a = 1.35 and the left-hand peak of these twin peaks is moving into the wall layer whilst 
the right-hand peak is decaying, showing the process by which the modes change up to the 
next higher wall layer mode. To complete the sequence, Figure 4e shows the eigenfunction 
for a = 2.5 and it turns out to be the eigenfunction of the fourth wall layer mode. 

Figure 5 shows solutions for a Mach number of 8 and we can see the first three adjustment 
layer modes. As the Mach number and wavenumber increase, the modes come very close to 
intersecting (although they do not as we shall show in Section 7) and this means we have 
to take prohibitively small steps in a in order to follow a particular mode numerically. As a 
result, some wall layer type portions of modes have been omitted from Figure 5. 


3. The Wall Layer Modes for M ^ 1 

Taking Equation (2.16) and using the scalings 

a 


a = 


M 2 


+ 


and 


13 + 


(3.1) 


(3.2) 
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allows us to look at the wall layer modes for M ^ 1 , for which 


V" + V' 
4 /" 2 


4/7" 


L(W' 2 )J 

1 


+ K 




(1 


'■zjst) J + jh - !)| (f (i - n + rr 


S' 

= 0 . 


(3.3) 


When a 1 we can get an approximation to the wall layer modes by the application of the 
WKB method. There are three regions to consider as we have a turning point associated 
with the WKB expansion. Close to the wall the transformation 


V{ V ) = S( V )(1 - f% 


(3.4) 


allows us to write (3.3) as 

s"( i - n - 2S(f f + r 2 + m + 5(1 - n [(- \a\i - i) 2 (i - f 2 ? 


r 


4 / 


7/2 


*_(r 

f \r 


r (1 -f 2 )) ■ 2 

which for a >> 1 gives the WKB approximation 


+ 5 ( 7-1 + 


= 0 , 


S" + a 2 S 
and we write this as 


1(7-1) (f“, 


-j(7 - i) 2 (i - n + [jo - n + f s 


S" + a J SIJI r/.fi) = 0 , 


(3.5) 


0, (3.6) 


(3.7) 


where 

H = -\{i - 1) 2 (1 - n + (j(! - f' 2 ) + f'f") > 

H > 0 for 0 < T] < T) t and r) t is the turning point. In order to satisfy the boundary condition 
(2.17a) we must have 5(0) = 0 since (1 — /') 7 ^ 0 for r) = 0. The solution of (3.7) which 
satisfies 5(0) =0 is 

S(rj) = -jr sin(a f (H ( 171 ))* drji) , (3.8) 

H 4 Jo 

where A is a constant and so V is given by 

V = sm(a r {HM)U m ). (3.9) 

H 4 Jo 

In the region of the turning point, the equation to be solved is the Airy equation 

V" — a 2 a{r] — rj t )V — 0 , (3.10) 
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where a is positive and given by 


a = — 


dH_ 

drj 




The required solution of (3.10) which decays exponentially as r/ ~> oo is 

V = Bit* Ai([d 2 a]z(r) - ? 7*)), 


(3.11) 


(3.12) 


where B is a constant and Ai denotes the Airy function (see Abramowitz and Stegun (1964)). 
If we denote the argument of the Airy function by t then since for t — 1 


Ai(t)~x *|f| * sin (J|f|§ + j) , 

we have that for t <C — 1 



Now letting r) — > rj t in (3.9) gives 


(3.13) 

(3.14) 


V = 


M i-/W) 

- v)* 



(i7(?/l))2j7/i). 


(3.15) 


Matching amplitudes and phases between (3.14) and (3.15) implies, since we allow B to be 
positive or negative in (3.12) and hence phases may be either in phase or i r radians out of 
phase, that 


A{1 - /'(^) 2 )a« 

1 

(J6 


= ±b, 


(3.16) 


and 



(3.17) 


where 

J_ = r(HM)idnu (3. i8) 

Jo 

and n is a large integer. 

In Figure 6 we have plotted the wall layer modes, with equation (3.17) for n running 
from one through seven and a Mach number of 8, and then superimposed them on Figure 5. 
We can see, even for a moderate Mach number and small values of n, that equation (3.17) 
gives a good approximation to the wall layer modes. 


4. The Adjustment Layer Modes for M 1 

For the adjustment layer modes, the appropriate equation below the turning point is 
(3.9). However, (3.12), the solution of the Airy equation (3.10), must now be replaced by 

V = Ett? Ai(t) + Fn? (4.1) 
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where E and F are constants since V will be large in the temperature adjustment layer away 
from the wall and so we are no longer seeking a solution which decays exponentially beyond 
the turning point. Thus for t <C — 1, since 

'2 


Bi(t) ~ 7r * |*| < cos Q|*|* + ^) , 


we have 


(E 2 + F 2 )* (2 3 3tt \ 

Vta ]ijl Sm \3 ^ 1 ~4 / 

where tan cj> = E/F. Matching (3.9) and (4.3) gives 

4(1 - f'{Vt) 2 )a« 


CT6 


= ±{E 2 + F 2 )* , 


and 

where n is again a large integer. 
When < > 1 we have 


_ n7T — -f (f> 


a - 


J. 


Ill 23 

Ai(t) = -7f _ 2r * exp (— — ta). 


and 


(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 


Bi(t) = IT h i exp (jt!) , 

so (4.1) gives for t 1 that 

i /2 A 

V & Ft * exp f -U ) , 

and this has to be matched with the solution above the turning point but below the tem- 
perature adjustment layer. In this region V is given by 

co(i -n 

exp^a 

J n t 


V ~ C °^ — ny l exp(a f {-H(rj i))*</i/i), 


(4.9) 


where we have only retained the exponentially large term, multiplied by a constant Cq, in 
the WKB approximation. Matching between (4.8) and (4.9) gives 


- = ^(i -WMn 

Co <76 

The temperature adjustment layer has z = 0(1) where 

'iogr + A 


t) = b + r — 


and b is a constant and 


r = ^log (Ml). 


(4.10) 


(4.11) 


(4.12) 
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In this region therefore 



/'« 1- 

M 2 ’ 

x 00 


(4.13) 

where h is a constant so 






T « 1 + (7 

— l)6e 2 , 


(4.14) 

and hence 






1)W * 

1(7- 

l)he z T 
M 2 ‘ 

OO 

(4.15) 


Therefore if we take z >• 1 and insert these expansions into (4.9) we get that the behavior 
of V as the temperature adjustment layer is entered from below is 


where 


V 


2cofc 


- 1)2 


e2* exp 


(aJ+ - 


l 


a (7 — 1 )he z | 

r Mi r 


(4.16) 


- 1: [fr - » ! <1 (f (1 - n + /'/")] *. 


2 M 2 f3 2 ' f v " J ' ' J J /I (4T7) 

In the temperature adjustment layer we make the transformation (4.11) in (2.16) and 
after scaling a and j3 by writing 

a = sr + . . . , 


0 = £rj + ..., 

we get, retaining the leading order terms for M ^ > 1, that 
1 — (7 — 1 )he 


<PV dV 

dz 2 dz 


1 + (7 — 1 )he z 


- 1 


+ Va 


~2 


-(1 + (7 - 1 )he’) 2 + i ( (7 2 1)feCa 


From (4.16) we can see that asz-too the behavior of V will be a multiple of 


s 2 e 


where 

5 = (7-I )he z . 

Letting 2 — ► —00 in (4.19) we can easily see that 


(4.18a) 

(4.186) 

= 0. 

(4.19) 

(4.20) 

(4.21) 


V ~ e 52 , (4.22) 

since we require V to decay above the temperature adjustment layer. Using (4.20) and (4.22) 
as boundary conditions we solved (4.19) numerically. In Figure 7 we have plotted the first 
three adjustment layer modes given by (4.19) for a Mach number of 8 and superimposed 
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this onto Figure 5. We can see that (4.19) is giving a good approximation to the adjustment 
layer modes. 

The numerical results we gained in Section 2 seemed to suggest the existence of a critical 
Mach number, above which the adjustment layer modes come into play and below which 
they do not exist, and we shall discuss this further in Section (3. 


5. The Wall Layer Modes for M '<» ~ 0(1) 


The modes located near the wall for ~ 0(1) will be described by Equation (2.16). 
For large a we are able to obtain an approximation to them from the WKB solution in a 
similar way to the solution obtained in Section 3 in the hypersonic limit. 

If we write 

/? = /? + a>, (5.1) 

then there is no turning point in the WKB solution. However, the solution obtained does 
breakdown when 7 ] ~ 0(a” 1 ) as a consequence of /' ~ arj for rj < Cl. Thus, in this region we 
write 

= aT w rj, (5.2) 

where T w is the value of T(rj) at 7 / = 0. (When the flow is adiabatic, i.e., the wall is thermally 
insulating, from (2.7), T w = 1 + (7 — l)M^,/2.) Substituting (5.1) and (5.2) into (2.16) we 
find that V satisfies the following equation: 


~/ 3 +2 



V 


(5.3) 


This is the equation satisfied by V in the corresponding incompressible problem described 
by Denier, Hall, and Seddougui (1991). It is a form of Whittaker’s equation and the solution 
for V satisfying the boundary conditions 


V = 0 at ^ = 0, 


and 


V — * 0 as — * 00 , 


may be given in terms of Kummer functions (see Abramowitz and Stegun 1964). Hence, we 
find that the unstable eigenvalues for the present case of a compressible fluid will be the same 
as those for the incompressible problem described by Denier. Hall, and Seddougui (1991). 
These are 

J_ 1 

’V2’ 73 ’"' ' 
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Thus from (5.1) the growth rates for the wall layer modes are given by 

_ n Cl 

13 = (5-4) 

where n = 1,2,3, — We note that the first solution corresponds to the exact solution 
described by (2.18). Moreover, the eigenvalues are independent of M 0 0 and T w . 

The appropriate solution of (5.3) is 

V(rJ>) — — n,2, 20), (5.5) 

where M(a, b , z) is Rummer’s function and b x is a constant. M{ 1 — n, 2, 2ip) is a polynomial 
of degree (n — 1) in 2 ip. The behavior of V for large is 

V(iJ>)~b 1 (-l) n - 1 ~e-U n + ... • (5.6) 

n\ 

Then from (2.16), we find that the WKB solution in the region above the wall layer which 
matches with (5.6) as rj — » 0 is given by 

V(v) = c i(TM) (1 -’ l>/ VMrexp[-a fTMJVil, (5.7) 

Jo 

where c\ is a constant. Since 

a 2 

f'(v) ~ otrj - —rj 4 + ... , (5.8) 

for t? <C 1, where a: = 0.332 . . ., c\ is given by 

ci =a n Ti ; 3n - 1 )/ 2 a- n (-l)"- 1 ^-6 1 . (5.9) 

In Figure 8 we show the growth rate given by (5.4) as a function of a for M — 5 for 
the first nine wall modes. These solutions are superimposed on the corresponding numerical 
solutions from Section 2. We see that for large a the asymptotic results are very good 
approximations to the numerical solutions. 

In order to discuss the near-linking of these wall layer modes and the trapped layer modes 
described in the next section, we must consider the solution of the wall layer modes for large 
mode number. This is apparent from Figure 3 since the growth rates of the wall layer modes 
will not be close to those of the trapped layer modes for a 1 unless n is large. 

Thus for n 1 we write the wavenumber as 

a — a$n + a\ + . . . , (5.10) 


12 


where cto and oq are constants. Then from (5.4) the growth rate of the wall layer modes for 
large n is given by 

dn «i 


/3 = ^ + ^V + .... 

2 2 8 2ao?i 

From ( 5 . 5 ) the solution for V in the wall layer for large n is 


_ 3 _ 1 _ 3 1 / I I I 3tt\ 

V ( = 0 i 2 4 7 r 2 ?z 4 ^ 4 cos ( 8 2 ^ 2 n 2 — 4 ”)* 


(5.11) 


(5.12) 


As ^ ► 0 the solution ( 5 . 12 ) matches onto the solution in a thin layer next to the wall of 

thickness 0 (n -1 ) where V satisfies 

d?V 2V 


d< 2 + T = °’ 


(5.13) 


and ( = nxf>. To satisfy the boundary condition at the wall and to match with ( 5 . 12 ) as 
£ — » oo we find that 

V ~ ej(( — ( 2 + . . .) forC<l, 

and 

3 i ^ / i.i 3tt \ 

V ~ bi2 4 x 2 n £ 4 cos ( 8 2 ( 2 as £ — > oo, 

where is a constant. 

For large values of xp the solution given by (5.12) is valid only for xp < 2 n and matches 
onto the solution above the wall layer for large n when xp is large and xp < 2n. We find that 
the WKB solution in the region above the wall layer which matches with (5.12) as r; — ► 0 is 
given by 


T 2 


V{tj) = gi{-Q(r})) ^T(r})cos((na 0 +ax) [ [-Q)h + 


where g\ is a constant and 6\ = 37 t/ 4 . The function Q(r]) is defined by 

n( \ <t( w* 2 (T(*l)f"(*l) T\v)\ 

qm = mi)) - - {—^ — j . 


(5.14) 


(5.15) 


and the solution (5.14) is valid for 0 < r) < rj* where rj* is the first zero of Q and Q < 0 for 
rj < rj*. A more thorough discussion of the function Q is given in the next section. 

In the region close to r] = r/*, from (2.16) V satisfies 


V" — t(tj — rj*)ar 0 n 2 V = 0, 


(5.16) 


where r = > 0. Since we require V to decay as ij — ► oo the appropriate solution of 

(5.16) is 

V = k 0 n*Ai(r), (5.17) 
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where k 0 is a constant and r = (aon 2 T)s (?/ — rj*). From the behavior of the Airy function Ai 
given by (3.13) for r — > — oo we have from (5.17) that 


V ~ ko\r\~* cos ^ 1 

(5.18) 

as r -> — oo. Then matching the amplitudes and phases of (5.14) as rj - 

-> rj* with (5.18) gives 

T ~ * ( a o n2 T )^ T(r]*)gi = k 0 , 

(5.19a) 

and 


(nao + — - — (7i + I 2 ) — 0i = — , 

ZCl 0 4 

(5.196) 

where the integrals I\ and I 2 are defined by 


h= 1 (~Q)*dr}u 

Jo 

(5.20a) 

and 


rv* 

h = / T 2 (-Q)-ld Vl . 

Jo 

(5.206) 

The expression (5.19b) determines a\, i.e., 


Oi + j- na 0 h 

ai ~ h-iiih + iiY 

(5.21) 


6. The Trapped Layer Mode for M^ ~ 0(1) 


For Mqo ~ 0(1) this mode is concentrated away from the wall and is the equivalent of 
the temperature adjustment layer made for Moo ^ 1 described in Section 4. For a 1 this 
mode can be regarded as a virtually continuous function of a. In Section 7.2, we discuss the 
near-linking of this mode with the wall layer solution for large n. Anticipating this we define 
a by (5.10) and consider the solution of the trapped layer mode for n 1. To obtain the 
solution for large a from the following results we set ai = 0 and a 0 n = a. For large n we 
write (2.16) in the form 


V" - 2=V" - a 2 n 2 QV = 0, 


where 


QM = r 2 


i - 


l ( f" 


fl 2 \Tf 

Then this mode is concentrated in the region where 


2T 2 , 


( 6 . 1 ) 


(6.2) 




Tf 2 T 


\T}=T) 


(6.3) 
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and 0 is a maximum at rj = fj. Therefore, from (6.2) we have 

Q(fj) = Q'(rj) = 0 and Q"(fj) > 0. (6.4) 

Thus when a 0 — 20q the function Q defined by (6.2) is also that defined by (5.15). 

In Figure 9 we plot Q(rj)/T from (6.2) for M ^ = 5 for a range of values of 0. We choose 
to show Q/T 2 rather than Q since (T( rj)) 2 given by (2.7) is large for small values of rj and 
decays exponentially to the free stream value of 1 as r/ increases. We see from Figure 9 that 
Q(v) — oo as T) — * 0+. For certain values of 0 the function Q(r]) has three zeros. For a 
particular value of 0 the second and third zeros coincide in a turning point. This position of 
rj is the location of the trapped layer modes r ] = fj, for large a. The sequence of modes will 
be described by corrections to (6.3). 

However, these trapped modes will not exist for values of below a critical value. We 
can see this if we look at the plots of Q(r))/T for various values of ft for = 3 shown in 
Figure 10. It is apparent from this figure that Q(r]) only has one zero for any value of ft and 
no local minimum for M 0 c = 3. For the Chapman constant C — 1 and an adiabatic fluid the 
critical value of M c 0 is found to be M ^ = M c = 3.564. Thus, for < M c the function Q 
will not have a local minimum and there is no solution described by (6.3) and (6.4). This 
explains why the discontinuous modes do not appear in the computational results for small 
Mach numbers. We note that there exists a critical value of the free stream Mach number 
when C ^ 1 or Prandtl number 1 and also for the case of an isothermal flow where the 
temperature of the wall is maintained at a constant value. 

Thus, from (6.3) ft is a constant to first order for the trapped modes and moreover, the 
same constant for each mode. The correction to ft will describe the distinct modes. In order 
to determine the eigenfunction in this region it is necessary to determine the correction to 
ft. In the trapped layer 

i i 

V = a%n*{ri - fj), (6.5) 

where rj 0(1). Write 

0 — 0o + a o ln 1 0i + • • • , ( 6 . 6 ) 

and substitute (6.5) and (6.6) into (6.1). Equating coefficients of the largest terms, which 
are of 0 (<XqU 2 ), shows that 0q is indeed defined by (6.3) and (6.4). 

Figures 11 and 12 show the solutions of (6.3) and (6.4) for an adiabatic fluid and also for 
the case of an isothermal fluid where the temperature of the wall is maintained at a constant 
value. In this case for Chapman constant C — 1 and a Prandtl number of unity, the basic 
temperature is given by 

T(n) = 1 + - fin)) 2 ) + (T, - 1)(1 - /'(>!)), (6.7) 
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where T w is the nondimensional temperature of the wall. Note that if T w = 1 + (7 - l)M^/2 
then (6.7) reduces to the expression (2.7) for an adiabatic flow. 

Figure 11 shows the value of /) for the trapped modes as a function of Mach number 
for an adiabatic wall and an isothermal wall with T w = 0.5, 1, 2, 4, 8. We see that rj 
moves towards the free stream as increases. For large values of the analysis must 
be replaced with that described in Section 4 so rj will correspond to the position of the 
temperature adjustment layer given by (4.11) and the modes will be as described in Section 
4. This is consistent with (4.11) since we see that the position of the temperature adjustment 
layer increases as increases. 

Figure 12 shows the corresponding values of fl 0 for the trapped modes as a function of 
Mach number. We see for the adiabatic case and for the isothermal case with T w < 8 that the 
value of /3 0 decreases initially as increases but then increases as M m increases further. 

The 0(a o n) terms in (6.1) give the following equation for V: 


cPV 

drj 2 



= 0. 

Po 


(6.8) 


where Q"(fj) is defined for /3 = (3 0 - 


The solution of (6.8) satisfying the boundary conditions 


V — ► 0 as |?)| — ► 00, 

is an eigenvalue problem for We find that the solution of (6.8) and (6.9) is 


(6.9) 


1/(5) = dMKif), 


( 6 . 10 ) 


where d\ is a constant, 

rj=(2Q"m*V, 

A ■ 2*(T(7?)) 2 /3 1 

(<?"(»?)) Mo 1 

and U(\,r}) is the Parabolic cylinder function (see Abramowitz and Stegun 1964). 
U(X,fj) grows exponentially as rj — »■ — 00 unless 


( 6 . 11 ) 

( 6 . 12 ) 

Now 


1 _3 _5 
2 ’ _ 2 ’ ~ 2 ’ ‘ 


Then the behavior of 77 (A, 77) as rj — ► —00 when A is given by (6.13) is 


(6.13) 


U(\ , 7 ;) ~ (— l)- A -i|^r A -ie-’ fS/4 , 


(6.14) 


and 




(6.15) 
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as rj — » +oo. 

From (6.12) and (6.13) we can evaluate ft for a fixed Mach number for a particular 
trapped mode. Figure 13 shows the growth rate from (6.6) for M ^ = 5 for the first three 
trapped modes as a function of a. In order to see the comparison between the numerical 
solution of the discontinuous mode described in Section 2, in Figure 13 only the discontinuous 
solutions are shown for large a and 0 below a certain value. We see from Figure 13 that for 
the first trapped layer mode the asymptotic solution is a very good agreement even for 0(1) 
values of a. For successive modes we still have good agreement but for larger values of a. 

Below rj = fj, but above rj = rj‘ where Q(ij*) = 0, Q(tj) will be positive. The WKB 
solution of (2.16) in the region r?* < r/ < fj is 


V{rj) = e 0 {Q{rj)) ^T(r/)ex p 


a 0 n + ai - 
V POj 


Q 2 drj i + 


ft n 

0o h* 


T 2 Q- 


drj\ 


(6.16) 


where eo is a constant. Now as r] — ► fj 

rv i rn L 

/ Q*drj i~ / Q 2 diji — a 0 n 
J rj* Jr;* 


7/2 

4 : 


and 


ft /Vn- A 


T'Q-UTjt-X ln((2Q"(f,))-* 


./Vq-u, 

Po Jr)* 00 Jr) 

where A satisfies (6.12). Thus, from (6.16) as rj — ► fj for large n 

t2 / 


2«ao n ^C t ?”(^)) *7\j))eolll * i(2Q")*exp 


J + (°° n + “ f") + 74 


(6.17) 


(6.18) 


(6.19) 


where /3 and are integrals defined by 


h 

n 1 

= / Q*drj, 
Jt}* 

(6.20a) 

and 

rrj-(aon)~% 2 1 

' T 2 Q~Urj. 


/, = tj 

(6.206) 

Po J 

n* 


Then (6.19) matches with (6.10) as fj — ► 

— 00 if 


2^aln^(Q"(rj))~^T(fj)(2Q l, )*e 0 exp 

^3 + ^4 = di( — 1) A 2 . 

(6.21) 


Close to the position rj = rj* V satisfies (5.16). We do not require V to decay above this 
region so the solution for V is 


V = koir 2 Ai(r) + l 0 ir 2 Bi(r), 


( 6 . 22 ) 
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where ko and lo are constants. 

From (4.8) the behavior of V as r — ► oo is 

V ~ lor-* exp • 

The solution (6.16) must match with (6.23) as t] — ► rj* . Then lo is given by 

lo = e 0 T~* (aln 2 T)^T(rf*). 


(6.23) 


(6.24) 


From (6.22) as r — * — oo 

v ~ (&o + cos + \ - 02 ) ’ ( 6 - 25 ) 

where tan 9? = ko/lo . In the region below rj — rj* V is given by (5.14), i.e., the solution above 
the wall layer for n 1. Thus matching the solution (5.14) as rj — > rj* with (6.25) gives 
instead of (5.19) 

T-Ha 2 0 n 2 T)uT( V ') gi = ±(k 2 + ig)l, (6.26) 

and 

(nao + a\)h — ^ — (h + I 2 ) — ^1 = O 2 — — . (6.27) 

Z&0 4 


7. The Near Crossing of the Wall Layer and Adjustment Layer Modes 

We see from the numerical results of Section 2 that for values of the free stream Mach 
number above a critical 0(1) value the wall layer modes and the adjustment layer modes 
become very close. This occurs for 0(1) values of the wavenumber and continues to occur as 
the wavenumber increases. 

A similar near crossing occurs for the inflectional acoustic neutral modes and the vorticity 
mode associated with the inviscid instability of a Blasius boundary layer for a compressible 
fluid and was discussed by Cowley and Hall (1990) and Smith and Brown (1990). The latter 
authors confirmed the conjecture of the former that, in the hypersonic limit, the modes are 
separated by an exponentially small amount. 

To discuss the near crossing of the two different types of mode described in Sections 3 
and 4 for M <*, 1 and in Sections 5 and 6 for a > 1 we follow Smith and Brown (1990) 

and extend our WKB analysis of these previous sections. 

Since the near crossing occurs for 0(1) Mach numbers as well as for M <*, 1 it should 

be possible to discuss the near crossing for both situations. The analysis will be similar for 
each case and therefore we attempt to minimize any repetition. We will first discuss the 
situation for the hypersonic limit. 
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7.1. The Hypersonic Limit 


We again have that below the turning point V is given by (3.9). In the region of the 
turning point V is given by (4.1) and (4.3) - (4.5) still hold. However, we replace (4.8) by 


i / » 4- 

t 4 


ji’exp (_!*!) + fexp (f<8)] , 


(7.1) 


for t 1. We also retain the exponentially small term in (4.9) so that V in the region 
between the turning point and the temperature adjustment layer is given by 


V 


(i-n 


c 0 exp|ay (-//(??i)) 2 d?;i| + d 0 exp |-a y (-//(??i)) 2 efyxj . (7.2) 


We already have (4.10) and matching the exponentially small terms between (7.1) and (7.2) 
gives 

5 §- = 4(1 - f(%n (7-3) 

ZCto CT 6 

so that 

£■=-■ p- 4 ) 

2 do Q) 

We now need to augment (4.16), the behavior as the temperature adjustment layer is ap- 
proached from below, by a corresponding term in d 0 so that 


2h* 


M 0 0 ( 7 -l ) 2 


T e ' d 


■{' 


,_ T a(7-l)h 

Co exp < aJ + — — — ttt e 


r ml 


! | + d 0 exp | - 


- T . a (7 - l)h , 

aJ + + 


(7.5) 


and then match this with the behavior of V in the temperature adjustment layer as z -► 00 . 
If we assume (4.19) has a solution 


V = V, (7.6a) 

a = A, (7.66) 

$ = B, (7-6 c) 

and then perturb these equations by putting 

V = V + V, (7.7 a) 

a = A + A, (7.76) 

p = B + B, (7.7 c) 
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and insert (7.7) into (4.19) we get after linearising 
1 — (7 — 1 )he 


<py_ dv_ 

dz 2 + dz [\l + (7 — 1 )he 

= V 


- 1 


+ VA 2 


~(l + (7-l )he Z ) 2 + ^( ( - 3 ^ 


2AA(l + ( 7 - l)/ie : ) 2 + (^r - <7 - !)*< 


5 


which can be integrated to give 


with 


w(z) = 2AA 


/;( 


V = Vw(z ), 


11 + (7 _D he^jy, d \^ 


(7.8) 

(7.9) 


F 


(a*b__ aa 
+ l B 3 B 2 


(7 ~^ h Io 


[1 + (7 — 1 )he Zl ] 2 [ z i 


F 2 


r 

Jo 


T 


>V 


[1 + (7 — \)he z,z \ 2 


dz 2 dz\ (7.10) 


Ci f z [1 + (7 — l)he Zl ]‘ 


(7-1 )h 


l 


F 


dz\ -f- C 2 . 


Using (4.22) we can calculate the constant C\ that is necessary to ensure that V decays 
exponentially as 2 — > —oo. We can then use (4.20) to calculate the behavior of w(z) as 
z -h ► 00 and if we choose C 2 so that there is no additional contribution of order s?e~ as to V 
from V then we have that as z — * 00 , V goes as 


V-[( 7 -l )he z Y 


exp{ — > 1(7 — 1 )he z ) + 


'-ac 3 - f~ - “pM - 


(7.11) 


where 


and 


a 


r — oc 0 roc 

= / v dz - 

Jo Jo 


V 2 dz , 


c t 


=1 


e z V 2 


-dz 


[1 + (7 - 1 )he z ) 2 
Matching between (7.5) and (7.11) gives us that 


L 


e z V 


[1 + (7 - 1 )he z ] 2 


dz. 


(7.12) 

(7.13) 


£0 

d 0 


AC 3 + 


' AB A\( 1 -l)h 


B 3 B 2 , 


C 4 


exp{— 2aJ + }, 


(7.14) 
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so we can get Fj E from (7.4). Assuming F is small since we want to perturb the wall layer 
modes (as well as the adjustment layer modes) for which F is zero, we get 



(7.15) 


and inserting this into (4.5) gives after some rearranging 



Setting the two brackets equal to zero gives respectively the wall layer modes and the adjust- 
ment layer modes. We can see from (7.16) that for finite these modes do not intersect 
but the distance between them is exponentially small. This explains the difficulties we had 
in following a particular mode when we solved (2.16) - (2.17) numerically. A sketch is given 
in Figure 14 of a pair of these near crossings with the continuous lines representing the situ- 
ation at a finite Mach number. The dashed lines represent the picture given by (7.16) with 
the exponentially small right-hand side ignored, i.e., when the adjustment layer modes have 
become continuous in the limit — * oo. 


7.2. Mqq ~ 0(1) 

As discussed in Section 5, in order to describe the near-linking of the wall layer modes 
and the trapped layer modes for M* ~ 0(1) we must consider the situation when the mode 
number of the wall layer mode is large. Similarly to the case when M ^ 1 we extend the 

WKB analysis of Sections 5 and 6 for n > 1. 

The solution in the wall layer for V is given by (5.12) with the solution below the position 
rj = if given by (5.14). The solution close to r? = rj* is given by (6.22) and so (6.26) and 
(6.27) are still satisfied. 

Now as r — > oo, from (6.22) we replace (6.23) by 

V~r~* yexp + foexp . (7.17) 

In the region rj* < rj < fj we must retain the exponentially decaying term so instead of (6.16) 
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we have 


^0?) = (QM) * T(r)) |e 0 exp |na 0 + a x - + j- 


Pi r-Ffd^-i 


+h 0 exp 


— nao + &i 




(7.18) 

where h 0 is a constant. Matching this solution as rj — » tj* with (7.17) gives (6.24) as well as 

(7.19) 

Following (6.19), from (7.18) as rj — * fj for n large 

V ~ 2t4>» i (C3"«))-J7 i K)(e 0 |f|- 1 -i(2(3"«))iexp 


ko ho 

2/o Cq 


V 2 ,( , M , ' 

-J + I flora + a-i - — 1 I 3 + u 


+/i 0 |r/| A ^(2 Q"{fj)) * exp 


T 

4 


Pi' 


a o n + o,\ — — 1 I 3 — I 4 
\ PoJ 


(7.20) 

Now we make the transformation (6.5) and perturb the solution in the trapped layer for 
large n where a is given by (5.10). Thus we write 


and 


P — Po H {Pi + P*) + • • • , 

a on 


V = F 0 + V*, 


(7.21) 


(7.22) 


where /3* and V* are small perturbations and V 0 satisfies (6.8). Then if we substitute (7.21) 
and (7.22) into (2.16) and linearize about the perturbed quantities we find that V* satisfies 


1 ‘ QPp-fv - Vh(T(fi)fv = 2 4-{T «)) 2 V„. 


(7.23) 


dr} 2 2 ' & '' v,// ' Po 

We wish to determine the behavior of V = V 0 + V* as 77 — ► —00. We want V*{fj) to 
decay as fj —1 00 so that the disturbance does not propagate outside the boundary layer. 
The solution of (7.23) is 

V*(r ) ) = tfO?)^), (7.24) 

where fj ) is given by 

* M) = f ( f mw rww** + D > 1; AAr + (7 - 25) 
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Using the behavior of V 0 as fj -> ±oo given by (6.15) and (6.14) we can determine the values 
of the constants D x and D 2 . Then from (6.15), since V* must decay for large rj we find 

A = /°>(M)) 2 <Si- (7-26) 

Po •'O 

For rj — > — oo we choose X) 2 to cancel the constant terms. Thus, from (6.14) we find that 

^(rj) ~ D 3 \rj\ 2X exp(rJ 2 /2), (7.27) 


as r] — > — oo where 

D. 3 = -(2 Q"{rj))-i 
Thus from (7.26), (7.28) is 

where 




D , 


A = ^(20"«))->TO)) 2 / 5i 


(7.28) 


(7.29a) 


too 

h= (7.296) 

JO 

Thus, from (7.22), (7.24), (7.27), and (6.14), the behavior of V as rj — > -oo is given by 

V(t}) ~ di(-l) _A_ 2[|^|- A -2 exp(-^ 2 / 4 ) + ^>3 1?7| A_ ^ exp(^ 2 / 4)] . (7.30) 


The WKB solution below the trapped layer must match with (7.30) as rj — ► — oo. Thus 
matching (7.20) and (7.30) gives 


£o 
h o 


(2 <rm 


D 3 


exp[— 2(na 0 + a x - ^-)h ~ 2 / 4 ]. 

Po 


(7.31) 


Thus, from (7.19) we have l 0 /k 0 . The solution for V described above is also a perturbation 
of the wall layer solution so we must have from (6.22) that lo is small. Then we have 

< 7 - 32 ) 

Z /Co 

Substituting this in (6.27) gives 

a.\{I\ — — — (/i + I 2 )) + naoh — ~ j ~T~- (7.33) 

ZCLq 4 rC o 


Using (7.19) and (7.31) we can rearrange (7.33) to give 

6 X + j - na 0 h (2 Q"(rj))~ 2 exp[-2(na 0 + a x - &)I 3 - 2 1 4 ] 
h ~ 2^(7i + h) 2 D 3 {h ~ j^(7 x + I 2 )) 


(7.34) 
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Now from (5.11) 


»-t 


a 1 


82 «q n 

Substituting the expression for a\ from (7.34) into (7.35) and rearranging gives 


°o 




1 /01 + 7r/4 — na 0 Ii \ 

ti4^\h-£;(h + h)) 




Po 


A 

a 0 n 


8foln(T«)p/ 5 (/,- 3 l ; (/ 1 + / 2 )) 


exp 


01 


-2 I na 0 + ai - — I / 3 - 2/ 4 


(7.35) 


(7.36) 


If we set the factors on the left-hand side of (7.36) to zero in turn we see that they 
describe the growth rates of the wall layer modes and the trapped layer modes respectively. 
Thus from (7.36) we see that for large n the wall layer modes and the trapped layer modes 
are separated by an exponentially small amount. This explains the apparent crossings of 
the growth rate curves obtained from the numerical results of Section 2. Thus the curve of 
growth rate as a function of a where the modes are very close is also described by Figure 14 
as for the case when 1. 


8. Discussion 

The main result of the present paper is that we are able to investigate the linear growth 
rate of Gortler vortices for 0(1) wavenumbers. This was achieved by considering the inviscid 
limit of large Gortler number for vortices having O(G^) growth rates. In this limit the 
growth of Gortler vortices is governed by parallel flow effects. Since for incompressible flows, 
or compressible flows with 0(1) Mach numbers, the spatial growth of viscous Gortler vortices 
is governed by non-parallel effects previous investigators were forced to consider the limit 
of large wavenumber, where as a result of boundary layer growth, non-parallel effects are 
unimportant (see Hall 1982a and Hall and Malik 1987). However, Hall and Fu (1989) show 
that in the hypersonic limit the growth of Gortler vortices with wavelength 0((21og M^,) - ^) 
is governed by a parallel flow theory. 

We have shown that in the inviscid limit, in a compressible boundary layer over an infinite 
cylinder, there are two types of growth rate modes possible for Gortler vortices. 

In the hypersonic limit we have firstly the type of mode we call wall layer modes. These 
are present in the incompressible case and their growth rates continue to rise as the wavenum- 
ber a — * oo. Secondly, we have the adjustment layer modes which have their eigenfunctions 
concentrated in the temperature adjustment layer away from the wall. The growth rate of 
these modes tends to a constant value as a — > oo. 
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We have also considered the limit a -> oo for 0(1) Mach numbers to describe the wall 
layer modes and the trapped layer modes (equivalent to the temperature adjustment layer 
modes when > 1). The wall layer modes are confined to a thin layer of thickness 0(a ) 

and have growth rates proportional to a* for moderate mode number n and tending to a 
constant value for large n. These results are the extension for compressible flows of the 
incompressible case described in Section 5 of Denier, Hall, and Seddougui (1991), hereafter 
referred to as DHS. For the compressible case we have an additional mode concentrated in 
a layer of thickness 0(«-|) away from the wall. This is the so-called “trapped layer mode,” 
which has growth rates tending to a constant value as a — > oo. 

The asymptotic results described for > 1 and also for M «, ~ 0(1) are shown to 
agree very well with the numerical solutions described in Section 2. We find that as the 
Mach number is increased from zero the wall layer modes start to deform until at a critical 

Mach number we see the appearance of adjustment layer modes. 

These modes are not present in an incompressible fluid. However, the growth rate of these 
modes tends to a constant value as the wavenumber increases while that of the wall layer 
modes (also present in incompressible flows) continues to increase. Hence, we anticipate that 
the trapped layer modes will not be as important as the wall layer modes for large values of 

the wavenumber. 

The situation here is similar to the case of inviscid disturbances to compressible flow over 
a flat plate. In this instance the vorticity mode does not appear until M*, exceeds a value 
approximately given by 2.2. However, in contrast to the present problem, the vorticity mode 

has larger growth rates than the acoustic modes when M <*> > 1- 

From the analysis for a > 1 for Chapman constant C = 1 and unit Prandtl number 
for an adiabatic boundary condition on the basic temperature we find that this critical 
Mach number is 3.564. For Mach numbers of this range, the adjustment layer modes are 
discontinuous with each mode solution of (2.16), apart from the exact solution, contributing 
a part to the adjustment layer modes as it passes through the relevant positions. 

We showed in Section 7.1 that for > 1 the wall layer solutions and temperature 

adjustment layer solutions of (2.16) come within an exponentially small distance of one 
another and in the hypersonic limit the adjustment layer modes become continuous. In 
Section 7.2 we showed that this is also the situation in the limit a -*• oo for M «> ~ 0(1) when 
the mode number of the wall layer modes is large. 

The Chapman viscosity law has been assumed in the present analysis, and results de- 
scribed for unit Chapman constant and unit Prandtl number for simplicity. There is no 
great difficulty in obtaining results for more realistic values oi C and Prandtl number and it 
is expected that the effect of this on the asymptotics will be small. 
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However, in the hypersonic limit the Chapman viscosity law does not give a realistic 
description of the viscosity of the fluid. Thus, of interest would be the results when the 
more accurate relation Sutherland’s law is used to describe the temperature dependence 
of the fluid viscosity. The results of Fu, Hall, and Blackaby (1990) who investigated the 
Gortler instability in the hypersonic limit for Sutherland’s law showed that the results were 
significantly different from those for Chapman’s law. In particular the structure of the 
temperature adjustment layer is altered, now having 0(1) thickness. See also, Blackaby, 
Cowley, and Hall (1990) who investigated hypersonic flows over a flat plate using Sutherland’s 
law. Thus, it is to be expected that significant changes to the results presented in this paper 
would exist if Sutherland’s viscosity law was used in place of Chapman’s law. 

The solutions of the equations governing the instability of Gortler vortices in the inviscid 
limit do not predict a fastest growing mode since the growth rates of the wall layer modes 
tend to infinity when a -* oo. From Section 5, in the high wavenumber limit we see that the 
wall layer modes have growth rates of 0(G2o§). We can extend the work of DHS to describe 
the fastest growing mode in the limit of large Gortler number for a compressible flow with 
M oo ~ 0(1). First we consider the structure of a viscous mode close to the right-hand branch 
of the neutral curve with a = \*G* . This mode exists in a layer of thickness 0(<z“?) centered 
around some non-zero value of 77. Initially A* is 0(1) but we are interested in the limit A* — ► 0. 
We note that the problem where A* -* 00 corresponds to the work of Hall and Malik (1989). 
The analysis is very similar to that described by DHS for the corresponding problem for an 
incompressible problem so the details will not be repeated here. As for the inviscid modes 
the growth rates are 0 (Ga). The first approximation to the growth rate 0 is determined 
from the following condition 



(¥♦¥?) 


Kuuy (u 0 , Ru 2 Ty (u 0 

~¥~ \T + —) + t T + ?A 


= 0, (8.1) 


and 0 is a maximum at y = y c where (8.1) is satisfied. Here JI is the viscosity of the 
fluid and T* is the Prandtl number. In the incompressible limit (8.1) reduces to precisely 
the expression given by DHS for the corresponding incompressible problem. We consider 
solutions of (8.1) when the basic flow is given by the compressible Blasius flow described by 
(2.4 - 7) for unit Chapman constant and Prandtl number. Then the situation of A* = 0 is 
precisely (6.3), determining the growth rate of the trapped layer mode. Thus, for A* = 0 
there will not be a solution of (8.1) for values of below a certain value. However, for 
A* > 0 we find that solutions of (8.1) exist for all values of We are not concerned with 
the details of the solution of (8.1) here but further particulars may be obtained, on request, 
from the second author. We are solely concerned with the behavior of the solution of (8.1) 
for A* — > 0 for the basic flow described in Section 2. The numerical solutions of (8.1) show 
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that 7 ] c becomes small for A* — > 0 and closer inspection of (8.1) in this limit reveals that 
7 ] c ~ A* 4 and /3 ~ A 1 "" 2 . Thus, in the limit A* 0 the growth rates of the viscous modes are 
0(G2(aG'4)- 2 ). Thus, as a decreases the viscous right-hand branch modes become more 
unstable. 

The above discussion suggests that there exists an intermediate region where the viscous 
mode described above and the inviscid mode of Section 6 overlap. This will occur when 
their growth rates are the same size. We find that this is the case when a ~ G$ with the 
growth rates of 0(G*). The vortices in this case will be confined to a wall layer of depth of 
0(G“* ). We note that this situation is identical to that described by DHS for the fastest 
growing mode in an incompressible flow. It turns out that the effects of compressibility 
may be scaled out with the result that the eigenvalue problem determining the growth rate 
of the fastest growing mode is identical to that solved by DHS. Thus we do not present 
the system of equations here but point out that full details of the compressible problem 
may be obtained from the second author. Consequently, the results of DHS also describe 
the solution for a compressible fluid. We have that each unstable mode has its maximum 
growth rate occurring at a finite value of A*. DHS showed that the most unstable mode 
corresponds to A = 0.476, j3 — 0.312 where a = AG> and the growth rate of the 

fastest growing mode is given by G^K^ct^T w 5 $. The results described above show that 
for a compressible fluid, as well as an incompressible fluid, the most unstable linear Gortler 
vortex at high Gortler numbers is viscous with wavenumber of O(G^), rather than 0(G*), 
which is appropriate to the unstable modes close to the right-hand branch of the neutral 
curve. An important result is that the most unstable modes occur close to the wall. This 
suggests that significant coupling coefficients will be possible in the receptivity problem for 
the most unstable modes. This was shown to be the case for the incompressible problem by 
DHS and identical results for a compressible fluid may be inferred simply from the results of 
DHS. Of interest would be the behavior of the fastest growing mode in the hypersonic limit. 
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Figure 6. The first seven wall layer modes given by (3.17) for a Mach number of 8 
(dashed lines) superimposed on Figure 5. 
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Figure 9. The function <3 (j?)/[T(??)] 2 from (6.2) for a range of values of 0 for M x - 5. 
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Figure 11. The position ij of the trapped layer modes for a > 1 as a function of M x for an 
adiabatic fluid and also for an isothermal fluid with T w = 0.5, 1,2, 4, 8. 





Figure 12. The neutral growth rate /? 0 from (6.3) for the trapped layer modes as a function 
of Moo for an adiabatic fluid and also for an isothermal fluid with T w = 0.5, 1,2, 4, 8. 


2 


3 


a 

Figure 13. The growth rate (3 as a function of a for =5: — numerical solution of (2.16); 
asymptotic solution from (6.6). 
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Figure 14. Sketch of the near crossing of the wall layer and adjustment layer modes. The con- 
tinuous curves represent the finite Mach number situation and the dashed lines the situation 
in the asymptotic limit — > oc. 
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